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Abstract 

I adapt a recently introduced method for integrating over the unitary group (S. Aubert and C.S. 
Lam, J.Math.Phys. 44, 6112-6131 (2003)) to the orthogonal group. I derive explicit formulas for a 
number of one, two and three-vector integrals, as well as recursion formulas for more complicated 
cases. 
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I. INTRODUCTION 



Integrals over compact Lie groups arise in physics naturally through applications of ran- 



dom matrix theory 1 



which plays an important role in various fields of physics, ranging 



fr0 ,„ nuclea, theory H , quantum chaos Q and transport in .nesoscopic dev.ces Q to q oan- 
turn information theory |5j. A new application arose recently in the latter field from the 
study of statistics of interference in quantum algorithms . Quantum algorithms can always 
be represented as real matrices by spending one more qubit (i.e. doubling the size of the 
Hilbert space), which labels the real or imaginary part of the wave function 0,0], and the 
question was posed how this influences the amount of interference necessary in the quantum 
algorithm compared to a unitary representation . Examining the statistics of interference 
drawn from the Haar orthogonal ensemble, i.e. the ensemble of orthogonal matrices with a 
flat distribution according to the Haar measure of the orthogonal group, then leads naturally 
to consider integrals of monomials of matrix elements over the orthogonal group O(N). 

Formulas for one- vector integrals were obtained earlier by Ullah 2\ and by Mehta jlfl| . 
An n— vector integral is defined as an integral containing elements from only n columns (or 
from n rows) of the orthogonal matrix. Prosen and co-workers introduced an asymptotic 
method which works well for large N 11 1. The problem was solved in principal by Gorin [12], 
who derived recursion relations connecting n— vector integrals to (n — 1)— vector integrals. 
While the method is general, we will see below that explicit formulas can be obtained in a 
much easier way in various simple but important cases of integrals of low degree (and/or a 
small number of vectors). A number of interesting properties of inte gral s over the orthogonal 
(as well as the unitary and symplectic) group were also derived in [l3|, based on the use of 
Brauer algebras. 

Recently, Aubert and Lam introduced a very elegant method for integration over the 
unitary group U(N) based solely on the unitary constraint and the invariance of the Haar 
measure under unitary transformations ^jj. This method is very powerful, and allowed 
the authors to obtain explicit formulas for the lowest order integrals, including all one- 
vector integrals, all two vector integrals with up to three different matrix elements taken 
to arbitrary powers, and all two and three-vector integrals up to order 6 (the order means 
here the number of matrix elements in the monomial to be integrated). Recursion relations 
can be obtained to reduce higher order integrals to these basic integrals. In the present 
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paper I adapt the method of invariant integration to the orthogonal group, and derive the 
corresponding basic integrals and recursion formulae. 



II. INTEGRALS OVER THE ORTHOGONAL GROUP 

Consider orthogonal, real N x N matrices O with matrix elements Oij, 1 < i,j < N. We 
will be interested in integrals of the type 



-L-llJl...lpJ p 



J(dO)O iljl ...O ipjp , (1) 



where the order p of the integral is some positive integer, and (dO) denotes the Haar invariant 
measure of the orthogonal group O(N), normalized to J{dO) = 1. These integrals can be 
calculated based solely on 

• the orthogonality relation T = 1 = 00 T (where T denotes the transposed matrix), 
or explicitly 

N N 

OijO ik = 5j k = ^ OjiO ki , (2) 



i=i i=i 
where 5jk stands for the Kronecker-delta; 

• and the invariance of the Haar measure, meaning that for any function f(0) and an 
arbitrary (real) orthogonal (N x N) matrix V, 

J (dO)f(O) = J (dO)f(VO) = J (dO)f(OV) . (3) 

Here we will use for f(0) the monomial 

f(0) = f[O i)Jx = I j, (4) 

A=l 

and we have introduced the notation Ojj, where I = {i±, . . . ,i p }, J = . . . ,j p }- 
Correspondingly, J,. ( ,.. will be abbreviated as T hh ,, Ap j v = l(0 L j). 

III. RELATIONS FROM INVARIANCE 

A number of powerful relations are obtained from choosing different explicit orthogonal 
matrices, under which Z(Ou) must be invariant. 
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A. Global sign change 



The simplest example consist of V — — 1, i.e. = — 8^, clearly an orthogonal transfor- 
mation, with O' = VO = —O. Eq.(j2J) implies f(0) = (— l) p f{0). Thus, p needs to be even, 
otherwise the integral will be zero. Even p will therefore be assumed from now on. 

B. Local sign change 

Another orthogonal transformation is induced by a matrix V with matrix elements Vij = 
(—l) Si 8ij, where Sj G {0, 1}. It leads to new matrix elements 0'^ = (VO)ij = (— l) Si Ojj and 
f{0') = (— l)^=i S4 a/(0). As each s ix is arbitrary 6 {0, 1}, eq.(j2J) implies that each i\ has 
to appear an even number of times. The same reasoning applied to multiplication with V 
from the right lets us conclude that also each j\ has to appear an even number of times, 
otherwise the integral will be zero. 

C. Permutations 

Consider the permutation that exchanges two indices, i$ <-> jo- It is induced by a trans- 
formation V with matrix elements Vij = 8^ for % ^ io,jo, V io j = Vj oio = 1, V ioio = Vj j = 0, 
which is manifestly orthogonal. An arbitrary permutation can be obtained by concatenating 
exchanges of two indices. The corresponding matrices are multiplied, and since the product 
of two orthogonal matrices is again orthogonal, an arbitrary permutation P can always be 
represented as real orthogonal matrix V. Multiplication from the left permutes the left 
hand indices in O, f{VO) = Y\\Op^ x )j x1 and multiplication from the right permutes the 
right hand indices, f{OV) = Yl\Oi\P(j\)- Therefore the value of the indices is of no impor- 
tance, the only thing which counts is the multiplicity of all different indices. Thus, we may 
rewrite / and J as / = {(i^ 1 , (« 2 )^ 2 , . . . , (i t )^} and J = {(j^ 1 , (j 2 )^ 2 , . . . , &)"•}, with 
Si=iA*i — V — Y^ii=i v ii wriere p is an even integer according to section IIII Al and all [ii 
and z/j must be even ( section Till B|) . We might even drop the indices i\ and j\ all together 
(e.g. chose them once for all as i\ = X (A = 1, . . . , t), j\ — A (A = 1, . . . , s)), and just keep 
the multiplicities M. = {/ii, . . . , /x t }, Af = {z/ 1; . . . , u s }. 

This result suggests a graphical representation of 1(0 u), where all first indices in Ou 
are represented by a dot in a left column (one dot for each different index) , and all different 
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FIG. 1: Example of a graphical representation for J (dO)OnOf 2 O2i022- Each dot in the left 
column presents a left index of a matrix element, each dot in the right column a right index of a 
matrix element. Each line corresponds to a factor in the monomial, with its power written next to 
it. The power is equal to one if it is not written out explicitly. The number of lines entering each 
point must be even. 

second indices in Ou are represented by a dot in a right column. These dots are joined 
by lines, where each line represents a factor Oi x j x . If a factor appears to the power m>, 
we denote that power next to the line (see fig. Q). If all powers m\ are even, the integral 
is positive, and is called a "direct integral"; otherwise it is called "exchange integral" 14]. 
In contrast to the integrals encountered for the unitary group [14], there are no complex 
conjugate factors here. This simplifies things at first, as we have to deal with only one type 
of lines, but also complicates things as it will lead to larger freedom in moving lines around 
(see below). Therefore, several recursion relations and their solutions will turn out to be 
more complicated than in the unitary case. 



D. Transposition 

We have (dO) = (dO T ), and therefore 1(0 u) = 1(0 jj). In the graphical language this 
tells us that all diagrams are invariant under reflection at a central vertical line. 



E. Rotations 



Consider a rotation in the plane spanned by the basis vectors pertaining to the indices 
a and b, = 8^ for i ^ a,b, V aa = V bb = cos£, V ab = -V ba = (sin£). Multiplied 
from the right, this transforms matrix elements according to 0' ia = (cos£)Oj a — (sin £)Oi&, 
0' ib = (sin £)Oj a + (cos£)Oj&, O'^ = Oij otherwise. It is most instructive to consider the 
effect of this rotation for a simple example, the integral I((On) d ), d even, and a = 1, b = 2. 



Expanding (OJ^, we obtain 



mO'n) d ) = g (^(cosO d - e (-sinO e | (d0)(0 n ) d - e (0 12 ) e = l((0 n ) d ) 



(5) 



where the last equal sign is dictated by the invariance under rotations, so the left hand side 
must be independent of £. As the only independent form of cos£ and of sin^ are powers of 
cos 2 £ + sin 2 £, the left side can only contain even values of e, such that with e = 2e, d — 2d 

l((0' n ) d ) = £ (If) (cos0 2{J - g) (sin0 2e X((O n ) 2 ^ g) (O 12 ) 2g ) (6) 

z— n \ / 



e=0 



^M e (cose) 2(d - S) (sm0 2 % (7) 



e=0 



^ = ('o~V((Oii) 2(d - S) (Oi2) 2S )- (8) 



2e. 

Invariance under rotations requires that M e oc The proportionality constant can be 
fixed from e = 0, M = T((O n ) fl! ) = Fi(2d), and we are thus lead to M e = (f)I((O n ) d ). 
Together with (jSJ) this leads to a series of relationships, 

F 2 (2(d-e),2e) = J((O n ) 2(J - g) (0 12 ) 2g ) = -|-X((On) d ) = ^F x (2^ . (9) 

In the graphical language, we have rotated 2e lines away from the right dot to a new, empty 
dot, and we have found a relationship between this new integral and the old one, containing 
only one line of multiplicity d = 2d. The same reasoning can be applied to a dot which is part 
of a more complicated diagram, i.e. a dot in which lines from several other dots terminate. 
The total number of lines arriving at the dot will again be denoted by d, and must still 
be even. Invariance under rotations implies again the same form (J2J) after expanding the 
powers of all matrix elements which get transformed, and leads to the same eq.([9*|). where 
now Fx(2d) means the original diagram, and F 2 (2(d — e), 2e) the diagram in which 2e lines 
have been rotated away from one dot to another, empty dot. Note that exactly the same 
relation is obtained for integration over the hypersphere (instead of integration over the 
orthogonal group) [l^. This is not surprising, as © is a one-vector integral, such that 
only the normalization of each column is involved, but no orthogonality relation between 
different columns. More generally, any one- vector integral over the orthogonal group O(N) 
equals the corresponding integral over the hypersphere in iV dimensions 12 1. 



FIG. 2: Graphical representation of the fan integral, ea. (|ll|) . 



Eq.Q can be iterated to give the relations called "fan relations" in 



Q, 



(^(i 1 )- ■( d - t - 2 ) 

F t (2(d- dO, 2(d! - d 2 ), 2K_ 2 - dVO, 2^0) = ^(2*-^^ 777^7 . (10) 

V2di/ \2d 2 / ' • • • ' l2dt_i/ 

We recognize 2d t ~i = m t as the multiplicity of the last line (line t), 2(d t _ 2 — d t _i) = m t _i 
as the multiplicity of line t — 1, and so on, up to the multiplicity mi = 2(d — di) of the first 
line. The relationships between dj and rrij are easily inverted, and, when re-injected into 
lead to the final form of the "fan relation" for an integral represented by the diagram 
in figE| 



„ , (d/2)! m\\ ■ . . . ■ mA 
F t (m u . . . ,m t ) = 1 > ^ . (11) 



2 r ■ ■ ■ \ 2 r 

The integral for a single line with multiplicity d, i*i(d), will be calculated in the next section. 



IV. RELATIONS FROM ORTHOGONALITY 

The orthogonality relation (|2*|) is a handy tool to reduce the number of lines in a diagram. 
It is most instructive to consider as simple example the integral 

F 2 (2(m-l),2) = |(rfO)(0 lfc ) 2 ( m - 1 )(0 lz ) 2 k^l, (12) 

As we have seen, the value of the integral is independent of the index I, as long as I 7^ k. At 
the same time, if we sum over all possible values of I (i.e. 1 = 1,..., N), we can make use of 
©> as EjIi(Ci0 2 = I- Therefore, 

(A-l)F 2 (2(m-l),2) = I \dO){0 Ve )^ m - 1 \O ll ) 2 - I \dO){O lk ) 2m 

1=1 

= F a (2(m-1))-F a (2m). (13) 
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This is an equation independent of the fan relation ([lip. The latter leads in the present case 
to 

F 2 (2(m-1),2) = ^, (14) 

and combining this with (|13j) gives a recursion relation for i<\(2m), 

2m — 1 

F 1 (2m) = F x (2m - 2) , (15) 

u ; 2m + A-2 u y 1 

with the obvious solution 

, s (2m- l)!!(iV-2)!! , . 

f '< 2m '= (2, n + iV-2)H • < 16 > 

The normalization -Fi(O) = J {dO) = 1 was used, and can be retrieved from (fTp^l if we define 
(-!)" = !■ 

The orthogonality relation is useful if we have a dot in which only two lines end. Note 
that so far we have used again just the normalization of each line in the orthogonal matrix. 
If the two lines ending in dot I on the right were replaced by two lines originating from two 
different dots on the left, the sum over all values of I would give zero, regardless of the rest 
of the diagram. 

V. Z INTEGRALS 

As an application of the techniques introduced, we now consider all possible integrals 
with up to three different factors, two indices on the right, and two on the left (see fieEIl), 
baptized "Z-integrals" in [14], 

Z( mi ,m 2 ,m 3 ) = J (rfO)(O n ) mi (Oi 2 ) m2 (0 22 ) m3 . (17) 

All multiplicities mi, m 2 , must be even. We start from the integral I(3b) shown in figOl 
obtained by rotating two lines away from the first dot on the left to a new dot. The latter 
is arbitrary, and we can sum over it, avoiding the already taken indices 1 and 2, whereas in 
the full sum the extra line disappears, 

N 

^J(36) =(JV-2)J(36) = Z{m 1 ,m2,m 3 -2)-Z{m 1 ,m 2 ,m 3 )-Z{m 1 ,m 2 +2,m 3 -2) . (18) 

fe=3 

On the other hand, the fan relation (jll)) can be applied to the upper two lines of the diagram 
(3b), which gives X(36) = m3 1 _ 1 Z(m 1 , m 2 , m 3 ). Once we insert the latter result into (JTHJ), we 
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FIG. 3: Graphical representation of the Z-integral, ea. (|17JI (left) and the integral I(3b) used to 
establish a recursion relation for it (right). 



find the recursion 

Z(mi, m 2 , m 3 ) = — ^ l — (Z(m 1 ,m 2 , m 3 - 2) - Z(m 1 ,m 2 + 2, m 3 - 2)) , (19) 

iV + m 3 — 3 

which should be supplemented by the initial value 

Z(m m 0)-F(m m ^ _ fa + m* - 1)\\(N - 2)\\ (^)!m 1 !m 2 ! 
Z(m u m 2 ,0) - F 2 {m u m 2 ) - ^ - - - - _ ^ ^ + ^j^,^, (20) 

obtained from the fan relation If we apply the recursion relation again, we easily 

convince ourselves that its solution is of the form 

Z(,n u m,, m 3 ) = ™f ^ z(nh , ,„ 2 + 2j , 0)( _ 1)( . (21) 

(iV + m 3 -3)!! ^ V « / 

A closed form can be found for the sum in this equation, if we consider N even and odd 

separately. For N even, (j2Dj) can be re-expressed as 

Z (mi, m 2 , 0) - — — j j - ma j j - mi+m2+jV _ 2 j j , (22) 



such that 



Zfm m m 1 - fa^SLlM S 
Z(m u m 2 ,m 3 ) - - - — - — ^^yS 1 , 

" ;; 2 7 -lY ////;; 2\ (m 2 • 2/)! 



4 y ^ j / ( m2 | z) ! ( rfl i+ rfl 2+2i+w— 2 j 



This leads to the final result 

Z{m u m 2 ,m z ) 



i=0 

2 vl2 ^{m 2 - l)!!(mi + m 3 + N - 3)!! 
(mi + iV - 3)n( ™i+^+™ 3+ — )! 



mi ! m2 ! m3! (m^)! (HMd^=2)l 



(23) 



2mi+m2+m 3 ( mi _|_ j\T _ 2)!(m 3 + iV - 2)! ( mi+r,1 3 +JV ~ 2 )! ^i+m 2 +m 3 +Af-2 ^j 

( mi +m 3 + ^-2)!(JV-2)! ' ' 

(m)!(m)!(^)! 1 ) 



2 2 -^r(^)r(^)r(^)r(iv - i)r(|(jv + m 1 + m 3 - i)) 

vrr(i(iV + mi - l))r(|(iV + m 3 - l))r(±(iV + mi + m 2 + m 3 )) 



(25) 
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The latter form turns out to be valid also for odd N. 



VI. EXCHANGE INTEGRALS 

Exchange integrals contain at least one line with odd multiplicity, and can therefore be 
positive or negative. The total number of lines arriving in any dot must, of course, still be 
even, otherwise the integral vanishes. Consider the structurally simplest exchange integral 
depicted in figlH 

x(r,s,t, u ) = J(do)(o n y(o 21 y(o 22 y(o 12 y, (26) 

r, s, t, u G N. The case r = s = t = u= lis easily solved by summing over the index of an 
arbitrary point, say the upper right one. We obtain 

^(M,M) = - F ^ T F 2 (2,2) = - (JV _ 1) ; (jv + 2) . (27) 

For general r, s, t, u we derive a recursion relation using invariance under a rotation between 
indices 1 and 3 (where 3 is a new index), multiplied from the right. This gives 



X(r,s,t,u) = X'(r,s,t,u) = J (dO) ((cos On - (sin0O 13 ) r ((cos0O 21 - (sin 0O 23 ) s 0* 2 

= (cosO r+s M + (cosO r+s " 2 (sin0 2 M 1 + ... (2! 

where M = X(r, s, t, u) is evident from choosing £ = 0, and Mi is given by the diagrams in 
figEJwith corresponding values X(4a), X(46), and Z(4c), 

Mi = r(r ~ 1) X(4a) + rsl(Ab) + S ^ S ~^ I(Ac) . (29) 

According to section IlII El we must have Mi = ?^X(r, s,t,u). The order of the integrals 
X(4a), X(46), and X(4c) can again be reduced by summing over the index of the new point, 
which yields 

X(4a) = j^-^ (X(r - 2, s, t, u) - X{r, s, t, u) -X{r-2,s + 2, t, u)) (30) 
X(46) = (X(r ' S ' t,u) + X (r-hs-ht+l,u+ 1)) (31) 

X(4c) = (X{r, s - 2, t, u) - X(r, s - 2, t + 2, u) - X(r, s, t, u)) . (32) 
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FIG. 4: The exchange integral X(r,s,t,u) (left) and the integrals T(4a), I(4b), and X(4c) used to 
calculate it (second from left to right, respectively). The latter integrals are obtained by rotating 
two lines away from dot 1 on the right to a new dot k ^ 1,2; the index k is then summed over. 

If we insert these integrals into ([29)1 . we get the recursion relation 

X(r,s,t,u) = -, -—- : -. : ; r ( — 2rsX (r — 1 , s — 1 , t + 1 , u + 1) 

v ''' ; (r + s)(N -2) + r(r - 1) + 2rs + s(s - 1) V V ' ' ' 1 

+r(r - 1) (X(r - 2, s, t, u) - X(r -2,s + 2,t, u)) 

+s(s - 1) (X(r, s - 2, t, u) - X(r, s - 2, t + 2, u)) ) . (33) 

No closed solution of this recursion relation could be found, but the principle of its appli- 
cation is clear: Iterating (J33|) will reduce the multiplicities of the lines with power r or s, 
till a power zero is achieved, and then the integral is reduced to a Z-integral, whose value is 
known for all multiplicities of the remaining lines, see eq.(|23|). Note that the formula corre- 
sponding to for the unitary case is quite different, as one must rotate always one line 
corresponding to a non-conjugated matrix element and one line corresponding to a complex 
conjugated matrix element (l4|. 

VII. INTEGRALS OF ORDER 6 

In this final section I will give explicit formulas for all integrals of order p — 6, which are 
not of the fan-type or Z-type. They are shown in figure El and will be denoted accordingly 
I (5a), I (5b), T(5c), T(5d), T(5e), 1(5 f), T(5g). All of them are obtained by summing over 
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-• •- 




• •- 




FIG. 5: All the integrals of degree p = 6 which are not of the fan- type or Z-type considered 
previously. These integrals are labeled 2 (5a), 2 (56), 2(bc), 2(bd), 2 (be), T(bf), T(bg) (from left to 
right), respectively, and can be obtained by summing over the indices of the dots with an arrow. 
2 (be) and 2(bd) are of the exchange type considered in fig0J 



the index indicated by an arrow. We are lead to 

1 



J(5a) 
X(5b) 
J(5c) 
J(5e) 
Z(5/) 



N-l 
1 

N-2 
1 

~N - 1 
1 



TV — 1 
1 

>-2 
2 

>- 2 



(F 2 (2,2)-F 3 (2,2,2)) 
(Z(2,0,2) - 2T(5a)) = 
F 2 (2,4) = 
^3(2,2,2; 



7V + 3 



(TV - l)7V(7V + 2)(iV + 4) 
-2 + 7V(7V + 3) 
(iV -2)(N- 1)N(N + 2)(N + 4) 
3(iV-2)!! 



(JV- l)(JV + 4)!! 
(iV-2)!! 



J(5d) 



(JV-l)(JV + 4)!! 



(J(5a) + J(5e)) 
J(5e) 



(N -2)(N -1)N(N + A) 
2(N-2)\\ 



(N-2)(N- l)(iV + 4)!! 



(34) 
(35) 
(36) 
(37) 
(38) 
(39) 



VIII. SUMMARY 

As a summary, I have adapted the method of invariant integration introduced in ^| to 
the case of integration over the orthogonal group O(N). Explicit formulas were obtained for 
all one- vector integrals (these coincide with integrals over a hypersphere in N dimensions), 
for all two-vector integrals with up to three different matrix elements, as well as for all 
integrals up to order 6. In more complicated cases recursion relations were derived, in 
particular for all exchange integrals with two different indices on the left and two on the 
right. 
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